An explicitly solvable model of the 
spontaneous PT-symmetry breaking 



in 

O ! Vit Jakubsky 1 , Miloslav Znojil 2 

O 

Abstract 

We contemplate the pair of the purely imaginary delta-function po- 
tentials on a finite interval with Dirichlet boundary conditions. The 
two parameter model exhibits nicely the expected quantitative fea- 
tures of the unavoided level crossing and of a " phase-transition" com- 
plexification of the energies. Combining analytic and numerical tech- 
! niques we investigate strength- and position-dependence of its spee- 

ch ' trum. 



^ : 1 The Model 

o . 

Systems with point interactions play an important role in quantum mechan- 
ics since they keep their exact solvability while describing real systems p. 
Conception of point interactions in pseudo-Hermitian quantum mechanics 
appeared in various works 0, [H], |3] 
O -1 The aim of present article is to describe spectral properties of one parti- 

cle confined in infinite square well with two point interactions of imaginary 
strength. The corresponding time-less Schrodinger equation reads 

c3 ; Hip = (- A - i£6(x + a)+ i£5{x -a))if) = Etj). (1) 

where H acts on continuous functions from L 2 ((— 1, 1), dx) vanishing in ±1 
and whose derivatives have an imaginary step ^//(±or)— ^//(±ol) = ±i£i/j(±a) 3 . 
Considering energy E as a general complex number, following the standard 
procedure we arrive at the secular equation 

sinh 2u H (sinh(4wa — 2u) + sinh 2u — 2 sinh 2ua) = 0, (2) 

Au 2 

where u is denned as E = —u 2 . Invariance of the equation with respect to 
complex conjugation of its root was rather expected due to the P-pseudo 
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3 Indexes L, R denote direction in which we approach the point. 
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Hermiticity of the hamiltonian (W = PHP). Due to the lack of physically 
relevant interpretation for non-real eigenvalues E, we will be interested in 
real energies only. In that case, (J2J) acquires the following simplified form, 

sin2/c + — sin2fcasin 2 k(l — a) = 0, k — Im(u). (3) 

K 

Neither (j2J) nor (jSJ) can be solved without numerical methods. Before pre- 
sentation of numerical solution, we will focus on qualitative properties of (J3j) 
that can be derived analytically. 



2 Spectral lines - robust vs. fragile 

Performing a formal manipulation, the equation (jSJ) can be rewritten, 

2 = sin2fc k 2 

sm2kasm 2 k{l -a)' 1 ' 

In (JTJ), we required £ to be non- negative real number. It restricts the domain 
°f £ = since the right-hand side of (jU must be positive. 
Let us take the following set 

^=(^n j iu(^n /+ )U M ' (5) 

where 

fc G 7 ± <£> ± sin 2A; > 0, k G J± ± sin 2fca > 0, 6 M <S> sin 2fc = sin 2fca = 0. 

Assuming that G J7", we define £ to be equal to its limit in the points where 
sin/c(l — a) =0. Then the function £ = is continuous on the open 
interval H-^- If ^ G M, equation (j3J) is fulfilled for any real £. 

Let us discuss now the behavior of £ at the boundaries of the subintervals 
of J. The numerator sin 2 A; equals zero in one boundary point at least. It is 
due to a < 1 and £ = in this case. We will say that the boundary point is 
of the first type. On the other hand when sin 2ka vanishes, £ tends to infinity 
and the boundary is of the second type. 

Thus, we can observe three kinds of spectral lines. The first one is sup- 
ported by an interval from J with both boundaries of the first type. It 
reaches its maximum ^ max on the interval so it is bounded. The spectral 
lines of the second type tend to infinity and are supported by subintervals 
°f D I T with mixed types of boundaries. The third kind is unbounded 
as well and is supported by k G M. That are semi-lines perpendicular to 
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Figure 1: The function £ = for a = 2/3. All three types of spectral lines 
are present, £ = 5t, k = J7r. 

fc-axes. In terminology used in the energy levels of the first type are 
called fragile while the remaining ones are robust (see Fig. 1). 

The third class of spectral lines corresponds to wave functions having 
nodes at the interaction points and it appears only for rational a. On two 
intervals j± C f] I T that are connected by k <G M, spectral lines of the 
first and the third kind are stuck together. We speak about unavoided 
level crossing in this case. Due to simple form of the level crossing can 
be described quantitatively. Let us take a = ^ where P and Q are coprime 
integers. The level crossing appears if 2k = IQtt and l(Q — P) is odd number, 
I is integer. In that case, the crossing of spectral lines occurs for 
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The above analysis provided us a qualitative knowledge about behavior 
of k = k(£) implicitly given by (J2J). The fragile spectral lines cease to be real 
when £ exceeds a critical value £ max . Thus for any a, we can find a sequence 
of the critical couplings £, C rit,n where n is the excitation of the higher energy 
from the couple of merging energy levels. The sequences of critical couplings 
have their minimum £ min ~ 2.4931 for a ~ 0.3335. It is an interesting feature 
of our model that the sequence is not monotonic. We call it an overstepping 
complexification of energy levels (see Fig 2.). The numerical solution k = 
k(£) of (0) in the supercritical regime is in Fig. 3. 



3 Conclusion and Outlook 

The presented model of a particle captured in an infinite square well with 
two point PT-symmetric interactions (0) offers an interesting spectral behav- 
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Figure 2: Overstepping complexification of energy levels for a = yj, £ = 8t, 
k = jir. 
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Figure 3: Real and imaginary part of E = — w 2 (£) after crossing exceptional 
point for a = |, w(0) = 7r. 



ior. We studied a parametric dependence of energy levels on a position and 
strength of the point impurities. The analysis provided us a good qualitative 
control over various phenomena appearing in the spectrum. We described 
presence of fragile and robust spectral lines as well as appearance of un- 
avoided level crossing of energy levels. These analytical observations were 
supported by numerical results. In the numerically obtained data, we ob- 
served a new phenomenon of "overstepping" complexification where merging 
of energy levels is not sequent with respect to the excitation but proceeds at 
random. 

The article is a continuation of research started in j^j. In this sense, it is 
straightforward to find perturbation series for energy in the case of generic a. 
It has been shown that for a = 1/2, the effect of non-Hermitian interaction 
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was negligible for high excitations. This fact enables approximate separa- 
bility of positive-definite metric performed originally in |7j and consequent 
computation of physically relevant quantities like position expectation value 
and its time development. 
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